gusucode.com > 《MATLAB图像与视频处理实用案例详解》代码 > 《MATLAB图像与视频处理实用案例详解》代码/第 21 章 路面裂缝检测识别系统设计/IterProcess.m

    function [bw, th] = IterProcess(Img)
if ndims(Img) == 3
    I = rgb2gray(Img);
else
    I = Img;
end

T0 = (double(max(I(:))) + double(min(I(:))))/2;
flag = 1; 

while flag
    ind1 = I > T0; 
    ind2 = ~ind1;
    T1 = (mean(double(I(ind1))) + mean(double(I(ind2))))/2;
    flag = abs(T1-T0) > 0.5;
    T0 = T1; 
end
bw = ind1;
th = T1;